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Abstract 

The shahow water equations describing the motion of thin hquid film on the rotating cyhnder 
surface are obtained. These equations are the analog of the modified Boussinesq equations for 
shallow water and the Korteweg-de Vries equation. It is clear that for rotating cylinder the 
centrifugal force plays the role of the gravity. For construction the shallow water equations 
(amplitude equations) usual depth-averaged and multi-scale asymptotic expansion methods are 
used. Preliminary analysis shows that a thin film of an ideal incompressible fluid precesses 
around the axis of the cylinder with velocity which differs from the angular velocity of rotating 
cylinder. For the mathematical model of the liquid film motion the analytical solutions are 
obtained by the Tanh- Function method. To illustrate the integrability of the equations the 
Painleve analysis is used. The truncated expansion method and symbolic computation allows to 
present an auto-Backlund transformation. The results of analysis show that the exact solutions 
of the model correspond to the solitary waves of different types. 

Introduction 

Many scientific and industrial problems connect to the flow of thin liq- 
uid films. Thin film technology is used extensively in many applications 
including microelectronics, optics, magnetic, hard and corrosion resistant 
coatings, micro-mechanics, biotechnology, medicine, laser, etc. At larger 
scales the ascent of buoyant magma below solid rocks and the spreading of 
lava on volcanoes are further examples of geological problems ||21] . Progress 
in these areas depends upon the comprehension of fiuid fiow mechanisms. 
As a rule the behavior of thin liquid films of an incompressible ideal fiuid 
can be described by the shallow water equations. Classic shallow water 
equations are obtained by the depth-averaging of the Euler equations for 
an incompressible fiuid under the assumption of potential fiow (see, for 
instance P-SIEIIT]). In this case the gravity plays the significant role. 
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The main objective of this paper is to construct the shallow water equa- 
tions for thin liquid layer (film) coated the surface of the infinitely long 
cylinder rotating with a constant angular velocity. In this case the role 
of gravity plays the centrifugal force. The usual depth-averaged technique 
and the multiscale asymptotic expansions method allows to obtain the ana- 
logues of the Boussinesq shallow water model and the Korteweg-de Vries 
equation. The main difference between classic equation and obtained model 
is the presence of the liquid layer curvature since the liquid motion occurs 
on the surface of a cylinder of finite radius. Naturally, with the tendency 
of a radius of cylinder to infinity the presented model pass into the clas- 
sical model. Other difference from the classical model is the presence of 
the vortex flow with constant vorticity (so-called the vortex shallow water 
equations, see for instance IHE])- 

Other objective of this paper is the investigation of analytic solutions for 
obtained equations. The study of exact solutions of the nonlinear PDEs has 
become one of the most important topics in mathematical physics. In the 
past decades, various powerful methods like the Inverse scattering method, 
variable separation approach and Homogeneous balance method were used. 
But, in recent years, much research works has been concentrated on the 
Cole-Hopf transformation method, the Jacobi elliptic method, Adomian 
method, and the various extensions of the Tanh-function method P-fT^ 

The paper is organized as follows. In Sec. [1] we introduce the basic 
equations governing the thin liquid film flow on the surface of a rotating 
cylinder. The corresponding two-dimensional free boundary problem is 
presented in Sec. [21 In Sec. El we apply the perturbation technique and the 
suitable transformations to construct a system of hyperbolic equations. In 
Sec. [H we derive the KdV equation on the basis of the Boussinesq model 
using ordinary technic of the amplitude equation constructing (see, for 
instance [8]). In Sec. [5] we present explicit Painleve test for the model 
equation. In Sec. [6] we solve the model equation analytically by using two 
different methods and discuss the fundamental properties of the model. 

1 The Basic Equations 

We assume that a thin layer of the ideal incompressible fluid with the free 
boundary coats the surface of the infinitely long cylinder rotating with a 



3 



constant angular velocity (see Fig.d]). To describe the behavior of this layer 
the Euler equation and the continuity equation rewritten in cylindrical 
coordinates (r, 6) are used 

1 v'^ 

Ut + UUr + -VU0 - — = -pr, (1.1) 
1 UV 1 , 

vt + uvr + -vve^ = — pe, (1-2) 

/j^ rj^ ry 

{TU)r + Ve = Q, (1.3) 

D = {a<r < R{e,t), < ^ < 27r}. (1.4) 

Here D is the region filled by the liquid, u is the radial velocity, v is the 
angular velocity, p is the pressure, a is the radius of cylinder. 
The equation of the free liquid surface is 

T = R{e,t), (1.5) 

where R{9^ t) is an unknown function that defines the free liquid surface. 




Figure 1: The sketch of hquid film which coat the rotating cyhnder 

The impermeability condition, the kinematic and dynamic conditions 
on the free boundary have the following form 

w = 0, r = a, (1-6) 

Rt + -vRe-u = {), r = R(9,t), (1.7) 
r 

P = u{e,t), T = R{e,t), (1.8) 

where I\.{9^t) is the pressure on the free boundary. 
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Naturally, no-slip conditions at the boundary of the liquid-cylinder are 
absent because of an ideal fluid. 

The dimensional and dimensionless variables are connected by formulae 
(dimensional variables marked by asterisk) 

{u'' .v"") = Q^R^{u,v), = p^QIrIp (1.9) 

Here Q^, U% R^, R^J^, p^, are the characteristic angular velocity, the 
angular velocity of the fluid, the characteristic radius, the inner radius of 
cylinder, the liquid density, and the kinematic viscosity, respectively. 



2 Fluid Flow with a Constant Vortex 

In this section we construct so called the vortex shallow water equations 
(see, for instance, P,[2]). We assume that vortex of the fluid flow uj equals 
to 2VL 

1 1 

UJ = -{rv)r U0, Lu = 2Q = const. (2-1) 

Obviously, the case of the potential fluid flow corresponds to Q = 0. 
We introduce the stream function 

u = —ri/jg, V = i/jr- (2-2) 

To determine the stream function ip we have the following underdeter- 
mined problem 

1 d dib 1 
= -jrrjr + ^^^^ = (2-3) 

^ = 0, r = a (2.4) 

It is easy to show that the solution of the problem ( 12.31) . ( 12.41) can be 
written as (see, for instance, [7] and Appendix 1) 

^(r, e, t) = -n(r^ - + E (-1)' (2/+1)! (^-^^ 

where F{6^ t) is an arbitrary function. 

We denote the stream function at the free boundary as 

"^{e.t) = ij{R{e,t),e,t). (2.6) 
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Then the kinematic condition at the free boundary (11.71 ) can be written in 
the following form 

RRt + = 0. (2.7) 

Using the dynamic condition at the free boundary (11.81) and (11.21) . (11.31 ). 

we get 

V f V v'^\ 

irv)t + uirv)r + -{rv)e + Re\ut + uur + -ue = -He, (2.8) 

r \ r r 



r = R{e,t). 

The system of equations (O), O with (O), dZSl), (ESD is a closed 
system of equations for determining the functions R{6,t), F{9,t). 

3 Long- wave Approximation 

Using the ordinary technique of constructing multi-scale asymptotic ex- 
pansions we introduce a fast variable 77, r and change scale of the variables 

T = et, T] = £0, u = eU, (3.1) 

where e is the small parameter related to the thickness of the liquid layer. 
In this case, equations (12 . 71) . ( ITHj ) . (12 . 21 ) . ( I?75| ) . (12 . 61) can be rewritten as 

RRr + = 0, (3.2) 

2 

V V 

(rv)r + U (rv)r + -{rv)„ - Rr, h 

+ £^R^ (Ur + UUr + ^/7^) = -n^, r = R{r], r), (3.3) 
*(„,r) = -a{R' - + E (-1)^' (2-TTy, d^P^'l'^)' (3.4) 

°o 2.7+1 n2j+l 

- = i'„ = E (-!>' e^' j^-^ ^ Fi^, r), (3.5) 
rv = ri,, = nr' + ^ (-1)^ e'' ^ F{n, r), (3.6) 



z = In ( - V Z = ln ^ 



M/ \ ^ / 

For convenience we rewrite (13. 41 ) -( 1X71 ) omitting terms of 
than 0{£^) 

^{tj, r) = ^n{R' - a') + ZF- s'^F^, + 0{e'), 

^3 

z^ 

TV = r^r = nr^ + F- e^^F^^ + 0{£^), 

z^ 

{rv)r = Fr-£^-F,,,,r^O 
Ur = - ^F,r + 0{£^), 



Ur 



1 



z 



ry 2 ry 2 

Z 



F,, + 0{£' 



r 



F^^ + 0{£'), 



3! 



(3.7) 
order more 

(3.8) 

(3.9) 

(3.10) 

(3.11) 
(3.12) 

(3.13) 
(3.14) 

(3.15) 

(3.16) 

(3.17) 



Substituting (I3.8I )-( IXT7| ) into (13.31 ) we obtain 



,Z 



Z^ 



F - E — F \ ZF - F — F 



.z 



2! 



R 



3! 



2^R-e'-F^A + 



|(s7fl^ + F-.^|!F,; 
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Retaining terms of order 0(£^) we have 
F — f'^—F —2QZF -^r'^—F F +2Qf^— + 



\ / \ J 2! 2! 



Finally, substituting (l3.8p -( l3rT7l ) into (13.21 ) and regrouping terms we get 
Boussinesq equations for describing the vortex shallow water 

RRr + nRR,^ + {ZF)^ = e^ (^^F,,,^ , Z = In (^^^ , (3.18) 

Fr+ (-2nzF -^n^R^ + nF + ^^] = (3.19) 



V 




/ rj 

If we neglect terms of order 0(£^) (i.e., the dispersion terms) then we 
obtain system of hyperbolic conserve laws 

(IrA + (l-nR^ + ZF] =0, (3.20) 



2 /_ V2 



Fr+ (^-2nzF-^n^R^ + nF + ^^^ =0. (3.21) 

Note that the investigation of this system has the independent interest. 
In particular, we can construct a self-similar solutions, rarefaction and 
shock waves, and solve the Riemann problem P,[3] 
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4 The Korteweg-de Vries Equation 

Using the ordinary technique for constructing of amplitude equation (see, 
for instance [8]) we can simpUfy the Boussinesq equations (13.181) . (13.191 ). 
We consider solution of these equations in the neighborhood of some char- 
acteristics of the linearized hyperbolic equations. Note that it is possible 
to construct models of different complexity using additional assumptions 
about the magnitude of the free surface perturbation. 

From our point of view we consider the most interesting case when the 
unperturbed free surface is a cylinder of radius c> a and the deviation of 
the free surface is sufhciently small (see fig.Ej) 



Here fih{r]^ r) is the function characterizing the deviation from the surface 
r = c, yti is the parameter characterizing the magnitude of the deviation. 
We also assume that the function F = 0{y) 



In general case e and /i are independent parameters. However, for sim- 
plicity we can link these parameters with help relation 



It means that magnitude of the free liquid surface perturbation is consid- 
erably less than the averaging thickness of the liquid. 




(4.1) 






(4.3) 




R 



c{l + iih{r],r)) 



The liquid film 



The unperturbed free surface 



\ The perturbation of the surface 

cfih{r], r) 



Figure 2: The deviation from the unperturbed free surface r = c 
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We express the function Z in the following form 

Z = lnf-) =Zo + /i/i + 0(/i2), Zo = ln-. (4.4) 
\ a a 



Substituting ( l4.1l )-( IOI ) into the Boussinesq equations (13.181 ). (13.191 ) and 
keeping only terms of order 0(£^) and 0(/i) we get 

c^(l + idh)hr + nc^{l + idh)h^ + ((Zo + idh)w)^^ = ^e'^Z^w^^^, (4.5) 



Wr + (— 2Q(Zo + ^h)w)j^ — Q^c^(l + ^h)hj^ + Qw,^ + ^wwr^ = 



= £^ Qw^,7,7r + ^^^W^,y,;,y^ - (4.6) 

Omitting terms of order 0(£^) and 0{jj.) we obtain a system of the 
first-order linear PDEs (not necessarily hyperbolic!) 

hr + + -4w„ = 0, (4.7) 

u;^ - n'^c^hr^ + 2nZo)wr^ = 0. (4.8) 

To determine the characteristic directions A of this system we have the 
equation 

{n - A)(Q - 2nZQ - A) + Q^Zo = 0. (4.9) 
Obviously, the characteristic directions (characteristic velocities) are 

Ai = Q(l - Zo) + n^Zi-Zo, M = 0(1 - Zo) - n^Zl-Z^. (4.10) 

We restrict our attention to the case 

Zo > 1. (4.11) 

Then the system (14.71 ). (14.81 ) is a hyperbolic system and there are two 
characteristics 

T] — AiT = const, T] — X2T = const. (4-12) 

We seek a solution of (14.51 ) , (14.61 ) in the vicinity of the characteristic with 
characteristic velocity Ai > 0. Using the asymptotic multiscale expansions 
method we introduce new variables (the index 1 is omitted) 

a: = 7/ — Ar, T = fir (4-13) 
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and the differential operators 

— - — - A— — - — f4 14) 

dr ^ dT dx' df] dx' 

In general case the functions /i, w are functions of the old and new 
variables, i.e. h = /^(r/, x^T), w = w{r], r, x, T). However, it can be shown 
that for a closed system of equations it is sufficient that these functions 
depend only on the new variables. 

We seek the solution in the form 

h = h^{x,T) + fih\x,T) + . . . , 

w = w\x,T) + iJ.w\x,T) + ... (4.15) 



For convenience we divide equation (14.51 ) on c^(l + fih) and omit the 
terms of order 0{fi^) 



hr + ^Ihrj + —Wtj + ^(1 - Zo)hWrj + = -s'^ ^Wrjrjrj, (4.16) 



2 2 2 2 

Wr — 2QZqWj^ — 2VLii[hw)r^ — Q c — Q c fihh^ + Qwjj H — 2^^?? — 



= £^ Qw^,;,^r + ^^W^rjri^ . (4.17) 

Substituting (KW^ into ( 14A6| ). ( ICTTj ) and using ( I4A4| ) we obtain 
(^^ - (/^"^ + /ti/.') + ^{h' + ^/i'). + §(^" + /ti^'). + 

+^(1 - Zo)(/i« + /i/i^)(w;« + iiw'), + ^(w;« + /iw;i)(/iO + iih% = 

= i/i%(w;V/iw;i),,„ (4.18) 
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- 2nZo{w^ + iiw^)^ - 2Q/i((/i° + iih'^){w^ + idw^)):, - 

H — ^(w + iiw ){w +iiw)x = 

{^^§T ~ ^"^^ ^ ^ ^ A^w;^),^,^ . (4.19) 

Collecting terms of the same powers of /i, we have 

-Xhl + nhl + ^w;; = 0, (4.20) 



-Xhl + r^/i^ + ^w;^ = (4.21) 



-Xwl - 2nZowl + nw;^ - n^c^hl = 



For obtaining the derivatives /i^, we have the linear system (14.201 ). 
The matrix of this system is degenerated and has the form 



12 



The left eigenvector (£1,^2) of the matrix M or the eigenvector of the 
adjoint matrix is given by relations 

(^Zq - vt^zl-z^ £1 - n'^c^£2 = 0, 



£i = n^c\ £2 = VlZq - VL^ Zl - Zq. (4.23) 
The solvability condition of the system ( 14.211 ) has the form 



/i^i + /2^2 = 0. (4.24) 

Integrating (14.211 ) with respect to x we obtain (arbitrary function of 
time T is omitted) 

»o = ;V.», N = ^^ = -^. (4.25) 

Zq Zq 



Using (14.251 ) and the solvability condition (14.241 ) we have 



h'r + -^Nh"hl + ^NhX - .^Nhl„ I, + (4.26) 



+iv - 4nhX + ^Nh^hl - nWhl + ^Aft°„ - InhlA e, = o. 

Finally, using formulae (14.231 ) for £1, £2 and (14.251 ) for N we obtain the 
Korteweg-de Vries equation (index is omitted) 

hr + aohh:, + Poh^xx = 0, (4.27) 

where 

£i{2 - Zo)c-^N - 4nN£2 + N^c-H2 - n^c'£2 

= £;tw2 = ^'-'^^ 

_ 3Q{Zo - ^Z'^-Zo){2Zo - 1 - 2^Zl - Zq) 



2Zo(l -Zo + ^JZl - Zo) 
„ -Zlc-^'Nh + 3AiV£2 - ^N£2 



-n(Zo - - Z^){-Zl + 5Zo + (6Z0 - 2) yzg - Zq - 6Zg) 
12Zo(l - Zo + ^JZl - Zo) 
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Now we return to the original variables (see, (13.11 ). (14.31 ). (14.131 )) and 
introduce the notations 



e = e-xt, x = n{i-Zo + ^jz^ - Zo), 

H = iih, R = c{l + H). (4.30) 
In this case the KdV equation takes the form 

Ht + aoHHe + ^oHeee = 0. (4.31) 

This equation describes the perturbation of the free surface in a coordi- 
nate system moving with velocity A and it can be used to characterize the 
nonlinear behaviors of the film flow traveling waves. 

In the fixed coordinate system it is possible to write 

Ht + XHe + aoHHe + l3oHe9e = 0. (4.32) 

The result obtained means that the KdV equation is a rough model of 
liquid film motion. Note, that we perform surface curvature effect with 
help constant coefficients ao and (3o. 

Other words, the previous Boussinesq model describes more accurately 
the behavior of thin liquid film on the rotating cylinder surface. However, 
the rough model also allows to describe the behavior of the liquid film. In 
particular, with help of the rough model we can describe the precessing 
liquid film motion in the azimuthal direction. 



5 The Painleve Analysis 

In this section we present the algorithm for the well-known Painleve inte- 
grabifity test which may greatly aid the investigation of integra- 

bility and the search for exact solutions. The more general definition used 
by Painleve himself requires all solutions of the ODE to be single-valued 
around all movable singularities. A later version [20] allows testing of PDEs 
directly without reducing them to ODEs. 

5.1 The algorithm and implementation 

We assume a Laurent expansion for the solution Hi 

oo 

//,(x) = 0-^(x)^//,,fc(x)0^-(x), ^z,o(x)7^0 and a. G Z", (5.1) 

fc=0 
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where Hi^k{'^) is an analytic function in the neighborhood of 0(x) . 

The solution should be single-valued in the neighborhood of the non- 
characteristic, movable singular manifold 0(x), which can be viewed as the 
surface of the movable poles in the complex plane. The algorithm for the 
Painleve test is composed of the following four steps: 

Step 1: (determination of the dominant behavior). To determine the 
strictly negative integer ai and the function i/^ o(x) we substitute 

F,(x) = //,,o(x)0"^(x) (5.2) 

into the following systems of M polynomial differential equations 

A(H(x),H'(x),H"(x),...,H(-)(x)) = 0, (5.3) 

where the dependent variable H has M components Hi, the independent 
variable x has A^^ components Xj, and H''"^^(x) denotes the collection of 
mixed derivative terms of order m. 

In the resulting polynomial system the equating of every two possible 
lowest exponents of 0(x) in each equation gives a linear system for deter- 
mination of ai. 

If one or more exponents ai remain undetermined then we assign a 
strictly negative integer value to the free ai so that every equation in (15.31 ) 
has at least two different terms with equal lowest exponents. Once ai is 
known, we substitute (lO) into and solve for H^™)(x). 

Step 2: (determination of the resonances). For each ai and Hifi{yi) we 
calculate the integer r for which iifj^r(x) is an arbitrary function in (15.11 ). 
We substitute 

i^,(x) = F,o(x)r (x) + i^,,,(x)r+^(x) (5.4) 

into (15.31 ). keeping only the most singular terms in 0(x), and require the 
equating of the coefficients iJi,r(x) to zero. It is correspond to determina- 
tion of the roots of det Q = 0, where the M x M matrix Q satisfies 

Q.Hr = 0, Hr = {Hi^r: ^2,r, • • • , Hu^r)^ ■ (5.5) 

Step 3: (determination of the integration constants and checking of the 
compatibility conditions). To possess the Painleve property the arbitrari- 
ness of Hi^r{^) must be verified up to the highest resonance level, i.e. all 
compatibility conditions must be trivially satisfied. 
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To verify these conditions we substitute 



//,(x) = r'(x) ^ //,,fc(x)/(x). (5.6) 

fc=0 

into ( 15.31 ). where VMax is the highest positive integer resonance. 

For the system with the Painleve property the quantity of the arbitrary 
constants of integration at resonance levels must be coincide with the quan- 
tity of resonances at that level. Furthermore, all constants of integration 
i^i^r(x) at non-resonance levels must be clearly determined. 

Step 4: (generation of the truncated expansion). For each pair of 
(oij, Hq{x.)) we calculate the possible truncated expansion in the form 

F(x) = Hor- + ^10"'+' + • • • + ^ = 1, . . . , M, (5.7) 

where Hi^{l = 0, . . . , —ai — 1) can be determined with substituting (15.71 ) 
into (15.31 ) and equating coefficients of the identical power of 0. 

If Hi, (/ = 0, . . . , —ai — 1) cannot be determined, then the series (15.11 ) 
cannot be truncated at constant terms. 



5.2 The Painleve test of the KdV equation 

To determine the dominant behavior, we substitute (15.21 ) into the KdV 
equation and remove the exponents of (j){6,t). The removing of duplicates 
and non-dominant exponents, and considering all possible balances of two 
or more exponents leads to 

ai = -2. (5.8) 

Substituting H{9,t) = Ho{e, t)(j)-^{e, t) into and solving for Ho{e, t) 

we get 

Ho{e,t) = -i2p4{e,t). (5.9) 

Substituting //(^,t) = -12/302(^, t)(/)-2(^, t)//,(^, t)0'-2(^, t) into ( I02| ). 
keeping the most singular terms, and taking the coefficient of Hr{9,t) we 
get 

r=-l,4 and 6. (5.10) 

While we are only concentrated on the positive resonances. The value 
r = — 1 is so called the universal resonance and corresponds to the arbi- 
trariness of the manifold H{9, t). The constants of integration at level j are 
found with help of the substituting (15.61 ) into (14.321 ). where tmqx = 6, and 
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by the removing of the coefficients at (j>^{6,t), such H^iO.t) and HQ{6,t) 
are arbitrary functions of t since the conditions at resonance r = 4, and 
r = 6 are satisfied. 

To construct the Backlund transformation of equation (14.321 ) we trun- 
cate the Laurent series according to Step 4 at the constant level term 

H{e, t) = Ho{e, t)4)-\e, t) + Hi{e, t)4)-\e, t) + H2{e, t). (5.11) 

Therefore, we obtain an auto-Backlund transformation of equation (14.321 ) 
as follows 

where H2 is a solution of the KdV equation. 

We take the vacuum solution at iif2 = in equation (15.111 ) which leads 
to 2 

H = Hor' + = 12/3 ^ (ln0). (5.13) 

Using the above auto-Backlund transformation and choosing the different 
H{6,t) and (f){6,t), one can obtain various solutions (as in |2D]). We can 
derive the special solution inserting (15.131 ) into (14.321 ) and collecting the 
terms with the identical power of (j). We get a system of homogeneous 
PDEs for 0. 

Finally, we get the following solution 

= C2 - 4/3ci) + i sin '^^^""'^ + sin 7(^ - 4/5ci)^ + C3, 

7 = (5.14) 

where Ci,C2 and C3 are arbitrary constants. 

This result means that the solution of the thin liquid film model corre- 
sponded to the Backlund transformation that is a solitary wave solution 
as shown in Fig. |3l 

Fig. [H shows the perturbation of the free surface of a rotating cylinder 
which describes the behavior of the thin liquid film. 
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Figure 3: Solitary wave solution of the KdV model equation by using the auto-Backlund 
transformation 




2 4 i E 2 4 S S 10 



Figure 4: The perturbation of the free surface of a rotating cylinder by Backlund trans- 
formation 



6 Analytical solution of the KdV equation 

To determine the solutions for various cases studied in this paper analytical 
schemes are used. The results for the nonlinear PDE ( 14.321 ) are presented 
and compared with the previous works and experimental situation p!2KT4] . 



6.1 Periodic solution of the KdV equation 

We introduce a self-similar variable ^ and the following notations 
i = - Qo = oioV, A = cotto, ^0 = oio(3. 



(6.1) 
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In this case equation (14.321 ) has the following form 



- VH' + cqH' + HH' + PH'" = 0. (6.2) 
Integrating (16.21) twice we get 

3/3{Hf = {hi - H){b2 - H){h - H). (6.3) 
Here hk are the constants of integration connected by relation 

V -cq = . (6.4) 

The periodic solution of (16.31) takes the form 

//(0 = ^dn2(z;s) + 63, (6.5) 

where 



Here dn^(z; s) is the Jacobi elliptic function 



2K{s) K{s) ^ l + g2n ^(5) 



9 = exp |-7r-^^^| , s=Vl-s2, (6.7) 

where i^(s) is the complete elliptic integral of the first kind. 

The function dn^(z; s) has a period 2K{s) and the wavelength A of the 
periodic solution is given by 

A = 2(^^^ ^ sK{s). (6.8) 

In the case of m waves on a circle we have 

A = — . (6.9) 

m 

The precession velocity V of the periodic solution and the average value 
H are given by 

V = co+'^-^^2_f^ + b,, (6.10) 

6S^ 
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A 

J? = A-'|ffK)d? = ||g + 63, (6.11) 



where E{s) is the complete eUiptic integral of the second kind. 

The minimum and maximum values of the function H{^) are defined by 
the formulae 

Hmin = &2, Hmax = ^1- (6.12) 

The periodic solitary wave solution (16.51) is shown on Figs. El 



H 




:100 



Figure 5: The periodic solitary wave solution of the KdV equation 




Figure 6: The perturbation of the free surface of a rotating cylinder 
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6.2 The Tanh- function method 

We use the transformation 

where the wave variable 

^ = e- Qot. 

This transformation converts the nonlinear PDE ( 14.321 ) to the equivalent 
ODE dO). 

The main idea of the Tanh-function method is to use a new variable 
tanh that allows to represent all derivatives of a tanh by a tanh itself. 
Introducing a new independent variable (and parameter u) 

Y = tanh[z/^] (6.13) 

leads to the change of derivatives: 

I - y^)^, (6.14) 

z/V 1 - y2 -2y-_ + 



I = - Y^) ((6y^ - 2) A _ ey(i _ y^)^ + (i _ yy^) . 

We construct solutions with help of the following finite series expansion 



m 



H{0 = S{Y) = Y,a^Y^ (6.15) 

in which all real constants should be determined later. 

The positive integer parameter m is obtained by balancing the linear 
terms of the highest order in the equation with the highest order nonlinear 
terms [22]. The highest degree of the linear term d^H/d^^ is taken as 

deg(^]=m + j, J = 1,2,3,... (6.16) 

and the nonlinear term H^d^ H/d^^ is taken as 

deg(^//^^^ =(fc + l)m + j. A: = 0,1, 2, 3,... (6.17) 
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We give a brief description of the tanh method as follows: 

Step 1 We determine the parameter m by balancing the highest-order 
partial derivative term and the nonlinear term in equation (14.321 ). 

Step 2 Using the standard mathematical (symbolic) or numerical soft- 
ware we substitute (16.151) into equation (16.21 ). equating to zero the co- 
efficients of all power y* yields an over-determined system of nonlinear 
algebraic equation for a^, v. 

We consruct the solution of the ID KdV equation (14.311 ) by using the 
abovementioned method. We can deduce from (16.161 ) and (16.171 ) into (14.321 ) 
that m = 2. It means that the KdV equation (14.321) may have the following 
traveling wave solution 



H{0 = ao + aiY + aiY^' 



(6.18) 



Substituting (16.181) into (16.21) . yields an algebraic system for ai{i = 0^1,2). 
The solutions to this algebraic equations can be derived 



ao = V — cq-^ 8/3/i^, ai = 0, 



09 



12/3 II' 



(6.19) 



These solutions represent solitary waves, such that the wave profile is pe- 
riodic, which gives good agreement with the the Painleve analysis and it's 
Backlund transformation as shown in Figs. [7], [HI 





Figure 7: The perturbation of the free surface of a rotating cyhnder 
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Figure 8: Solitary wave solution of the KdV model equation by using the Tanh-function 
method 




Figure 9: The perturbation of the free surface of a rotating cylinder 

On the Fig. [7H2] a rich dynamical behavior characterized by appearance 
of nonHnear waves is iUustrated. The motion of the film in two spatial 
dimensions with a rotating cylinder and without surface tension exhibits 
a solitary waves that corresponds to the decomposition of the nonlinear 
waves in previous models |[I2HI1]. 



Conclusion 



In this work we have obtained three models, describing the behavior of 
a liquid film of an incompressible ideal fluid on the surface of a rotating 
cylinder. Boussinesq model (13.181) . (13.191) . the most accurate one, takes 
into account the mean curvature of the liquid layer. The equation of KdV 
(14.321 ) is a rough model, which allows producing records due to a curvature 
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of constant coefficients. Tfie system of quasilinear equations ( 13.181 ). ( 13.191 ) 
in the case of hyperbolicity also allows to obtain information about the be- 
havior of the liquid film. Despite the fact that the model ( 14.321 ) is rough, it 
nevertheless allows us to construct a solution corresponding to the preces- 
sion thin film along the azimuthal direction, determine the shape of the free 
surface (in the form of cnoidal and solitary waves) and the rate of preces- 
sion. We have applied directly the PDE Painleve test to the KdV equation. 
An auto-Backlund transformation is presented by using of the truncated 
Painleve expansion and symbolic computation. The Tanh-function, Jacobi 
elliptic function expansion methods are used to solve the KdV model. The 
solutions describe the perturbation of free surface. At the onset of this 
phase transition, the system reveals rich dynamic behavior characterized 
by the appearance of solitary waves. In addition, we find that the solu- 
tion which is obtained from the Painleve analysis behaves similar to the 
solution of the KdV equation with the Tanh-function method. 

Appendix 1. The derivation of (12.51) 

To determine the function tjj{x, z) we have the following problem 

+ ^zz = ^ £^ tl^ = ij{x,z,t), (A.l) 
V^(x,0,t) = 0, (A.2) 

where Q is the constant. 

We seek a solution in the form 

oo ^ 
^Pix, Z,t) = Y^ '0fc(x, Z,t) + -VL£^ (A.3) 
A;=0 

^fc(x,0,t) = 0, /c = 0,l,2,... (A.4) 
Substituting ( IA.3I ) in ( lA.ll) and collecting terms with identical powers 



of £ we obtain 

ipQ^zz{x,z,t) = (A. 5) 

^k,zz{x,z,t) = -'ipk^i.j,^{x,z,t), A: = l,2,... (A.6) 
Taking into account the boundary conditions ( IA.4I ) we obtain 

ipQ{x,z,t) = co{x,t) z, 

^fc(x, z, t) = Xl(-l)'^' (2l-l)\ ''k+l-^{x, t), k = l,2,..., (A.7) 

i=l ^ ' 
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where Ci{x,t) are arbitrary functions. 
Then function i/j has the form 



oo 



fc+1 



fc=0 



i=l 



2i - 1)! 



i-2 



1 
2 



oo fc+l oo oo 

By changing the order of summation ( = X] X ) obtain 

k=0 i^l i=l k=i-l 



00 



' ' ^ Z^v ; 27 + 1! 5x2.?^ 1^ jy ^ J 2 

J=0 ^ ^ fc=j 



Finally, we have 



oo 



2j+l n2j 1 

J=0 



(2j + l)! 5x2^' 



(A.8) 



where 



00 



(A.9) 



i=0 



Of course, the formula (I A. 81 ) is formal. At least, this formula is valid 
at e = 1. It is easy to verify with help of direct substitution of (I A. 81 ) into 

(EH), (D). 

By analogy, we can obtain a similar formulae for the polar coordinate 
system. 

Let we have the following problem 



r or or r^ od^ 

t/j{a,e,t) = 0. 

We introduce the changing of the variables 

r d d 
^ = ln-, r— = — . 
a or oE, 



(A.IO) 
(A.ll) 

(A.12) 



In this case the formulae (lA.lOp . (lA.lip take the form which coincides 
with (EH), (ID) 



(A.13) 
(A.14) 
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Obviously, the solution of the problem ( lA.lOp . (lA.lip is 



^(r, e,t) = J2 i-^y^'' -^^^^^ ^ Hie, t). 



(2j + l)! o»(92j 



(A.15) 



Appendix 2. Jacobi elliptic functions 

We present a set of formulae used in Sec. IGT , 

dn^(/it) = 1 — sn^ {fit), 



d 



d 



— dn (/it) = —2s sn(/it)/i- — sn (/it), 



(it 



(i/it 



sii(/^i) 



fit 



dz 



d 

dfit 



sn(/it) = V(l - sn2(/it))(l - s2 sn2(/it)). 



|^dn2(/it)| = V(l - dn2(/it)) dn2(/it)(s2 - 1 + dn2(/it)) 
Let 

u = adn^(/it) + c. 
Then m satisfies the equation 



a 



u'Y = {a ^ c — u){c — u){a ^ c — as''' — u). 



4/i2 

We introduce notations 



a + c=6i, a + c — as^ = &2, c = 63, 3/3 



a 



4/i2 



In this case equation (IB.7P becomes 

3^(w')' = (6i-w)(62-w)(63-w). 
Also, we have the relations 



26 61-62 2 ^1-^2 

— = 61 - 63, 6 = — - — , s = -. 

2 61 — 63 



:b.io^ 
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